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When a fluid jet strikes an inclined solid surface at normal incidence, gravity creates a 
flow pattern with a thick outer rim resembling a parabola and reminiscent of a hydraulic 
jump. There appears to be little theory or experiments describing simple aspects of this 
phenomenon, such as the maximum rise height of the fluid above the impact point, and 
its dependence on jet velocity and inclination angle. We address this with experiments, 
and present a simple theory based on horizontal hydraulic jumps which accounts for the 
rise height and its scaling, though without describing the shape of the parabolic envelope. 
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1. Introduction 

The description of the hydraulic jump arising from a jet striking a solid surface has been 
a rich source of fluid dynamical problems for decades. The first fairly complete description 



is usually attributed to Watson (1964), and his theory has been refined and improved 



by many authors (for example Bohr et al. (1993); Godwin (19931; Bohr et al. (1997); 



Brechet fc Neda (1999); Chang et al. (20011; Bush fc Aristoff (20031). Experimentally, 



there are several variants of the problem, such as where the jet strikes a horizontal plate 



gle ( 


Sparrow & Lovell|1980 


Rubel||1981 


Kate et aZ.||2007) 


(|Gradeck et aZ.||2006 


Kate et aZ.||2009 


)■ 



Another important case involves the impact of a jet on a plane inclined with respect 
to the horizontal, so that there is now a gravitational force tangential to the plane. A 
profile resembling a parabola is then observed, with a maximum 'rise distance' along 
the plane (figure [T]). This rise distance can be associated with the hydraulic jump in the 
purely horizontal case, although there are differences. A full description of this problem is 
challenging, partly because of the breaking of axial symmetry, but also due to the neces- 
sarily complex flow in the thick 'rim' bounding the inner, shallow region. For that reason, 
most existing models are based on inviscid fluid dynamics (or the opposite extreme of 



Stokes flow, as in Lister ( 1992 )). Rienstra ( 1996 1 first described this situation in terms of 



the ballistic motion of fluid particles, leading to parabolic trajectories with a parabolic 
envelope or rim. Such a ballistic model predicts a rise distance of order U 2 /(2gsm9), 
where U is the jet velocity and gsm9 the component of gravity along the plane. In this 
model the particle trajectories are characteristics of a geodesic equation modified by grav- 
ity (Thiffeault & Kamhawi 2008), and these characteristics cross. Edwards et al. (2008) 



used a 'delta-shock' model to resolve the crossing of characteristics, which lowers the rise 
distance to 5/9 of the ballistic value whilst maintaining a near-parabolic outer envelope. 




However, in their inviscid model the rise distance still scales in the same way as for the 
ballistic theory. In their paper Edwards et al. ( 2008 ) performed one experiment (at 90° 
inclination angle) and found a significantly lower value of the rise distance than their 
theory predicted; they expressed hope that faster flow rates might capture the inviscid 
regime, but this has still to be verified. 

In the present paper we account for the discrepancy using a viscous theory, motivated 
by observations made in a simple experiment. The theory is a straightforward modifica- 
tion of Bohr et al. (1993), where the authors matched inner and outer solutions to viscous 
shallow-water equations to predict the radius of a circular hydraulic jump. To simplify 
the treatment, we include a component of the gravitational force pointing towards the jet 
while maintaining the assumption of axial symmetry. The resulting equation correctly 
captures the dependence of the rise distance on both the jet velocity and inclination 
angle. The model, however, is not sufficient to address other features of the hydraulic 
jump, such as the shape of the envelope, or whether the jump closes or has an open (i.e., 
parabolic-like) shape (see Lebon et al. (2008)). The model does not agree as well with 
experiments at slower rates of flow, as in Bohr et al. (19931, or at angles of inclination 
larger than about 60°. 



2. Experimental setup and results 

We performed a simple experimental study on an inclined hydraulic jump to measure 
the dependence of the rise distance on flow rate and inclination angle. A large plexiglas 
sheet is held clamped over a sink at a constant angle, so that the water runs directly off 
at the edge. This sheet is painted white on one side and marked with a regular 10 cm 
grid for length calibration. A small pump is used to supply fresh water at a constant 
flow rate Q, connected to a straight glass tube by flexible plastic tubing. The glass tube 
defines the nozzle of the jet, which has an inner diameter d = 0.56 cm. It is fastened to a 
metal rod that can be adjusted so that the jet strikes the sheet at normal incidence for 
each value of Q and inclination angle, while the nozzle itself is kept at a 3.5 cm distance 
from the sheet. The flow rate was measured after each run by the time taken to fill a one 
liter beaker for each Q. Based on this and the nozzle geometry we have an exit velocity 
U ~ 1-4 m/s, with Reynolds number Re ~ 6 x 10 3 -3 x 10 4 . 
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Figure 2. Maximum rise distance i?j ump , measured upslope from the center of the jet t o the 
lower edge of the bounding rim. The solid lines come from the theoretical prediction (3.141. The 
theory works well for moderate angles (not near 90°) and faster jet speeds. 



We inclined the plexiglas at five different angles 9, measured with respect to the hor- 
izontal: 9 = 6°, 15°, 28°, 45°, and 90°. At each angle, we varied the flow rate Q ~ 20- 
150 cm 3 s" 1 . Finally, from photographs we measured the rise distance i?j U mp, which we 
define as the distance from the center of the jet to the position of the zenith of the hy- 
draulic jump (i.e., the lower part of the rim or envelope at its point of highest rise). Note 
that at higher Q and 8 the flow becomes more unsteady, and the rise distance has larger 
error bars (in those cases we average the rise distance over time). Figure [2] summarizes 
the results for the measured i?j U mpi which increases with velocity and decreases with 9, 
as one would expect. We now present a simple theory to explain this dependence. 



3. A simple model 

To capture the maximum rise distance of the fluid along the sloping surface, we will 
model the flow as axisymmetric with a radial force, where the radial force arises from 
gravity along the slope angle. The radial model is suggested by the radial surface waves 
visible in figure [T] which appear roughly circular despite the inclination. This simplified 
model will prove sufficient to capture the scaling of i?j U m P with jet velocity and inclination 
angle. To satisfy mass conservation, one can imagine truncating the surface just past the 
hydraulic jump, so that fluid can spill out. This model resembles a truncated inverted 
cone, with the jet hitting the apex, though without the geometrical factors associated 
with a cone. 

Following Bohr et al. (1993), we start with the steady, incompressible Navier-Stokes 
and mass conservation equations in axisymmetric cylindrical coordinates, in the boundary 
layer approximation: 

du du dh . d 2 u . . 

u- — \- w-—- = — gcosti- qsmv + v—. r, (6.1a) 

Or Oz dr az z 

du u dw .„ „ , . 

^ + - + ^=0, 3.16 
or r Oz 

where u(r,z) and w(r,z) are respectively the velocity components tangent and perpen- 
dicular to the solid surface, g is gravitational acceleration, h is the fluid thickness, and 
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v is the kinematic viscosity. The coordinate r is tangent to the solid surface, and z is 
perpendicular to it. The boundary conditions at the bottom and top of the fluid are 



du 
dz 



0, w = 0, 
dh 



0, w 



dr 



at z 
at z 



0; 

h(r). 



Integrating equation (3.16) gives the mass conservation equation 

eh 

u(r, z)dz — q 



(3.2a) 
(3.26) 

(3.3) 



where q = Q/2ir, with Q the flow rate of the jet. In addition, we must specify the 
velocity u(ro, z) = uq{z) and height h(ro) — ho at a radius ro larger than the jet radius, 
since the boundary layer equations are not valid directly under the jet. When 6 = 0. 
equations ( |3.1[ ) reduce to those of Bohr et al. (19931. 

We use the standard hydraulic jump scalings to define dimensionless 'tilde' variables: 



it 
w 
r 



cm, 
/3w, 
Tf, 



a 



= Si 



{c- l/2 c y & )q x l s v x '\g C os6fl\ 

( Cl /2 C 2 - 3/8 )Q 5 /8 l ,-3/8 (5C0s ^-l/8 ) 

qVW^gcosd)- 1 / 4 , 



S = icT) 



(3.4a) 
(3.46) 
(3.4c) 
(3.4d) 



except that we included the cos 9 dependence in the scalings. (We will discuss the dimen- 
sionless constants c\ and ci below.) We immediately drop the tildes. 



In dimensionless variables, the mass conservation equation (3.3) becomes 

rh 

u(r, z)dz = 1. 



(3.5) 



Upon averaging in the z direction, equation (3.1a) becomes after integration by parts 



2u 



du + 1 
dr r 



= —ci G 



dh\ 
dr) 



ci 1 du 

C2 h dz 



(3.6) 



2=0 



where the ^-average of a function F(r, z) is F(r) = (1/h) / ( F(r,z)dz, and we defined 



1/2 -5/8n 3/8 -5/8 



(g cos ( 



a/8 



tan/ 



We now assume the separable form 

u(r,z) = v(r)f'(z/h(r)), 



(3.7) 



(3.8) 



where v = u is the averaged profile, / is a given function that describes the vertical 
structure of the thin layer, with /(0) = f (0) = /"(l) = 0, /(l) = 1. With this form for u, 
the mass conservation integral (3.5) becomes simply vhr = 1, which gives a relationship 
between v and h. This allows us to derive the two relations 



2 a 



du 
dr 



ti 



C\ v 



dv 

dr ' 



du 



C2 



2 = 



with 



ci 



and 



C2 



/"(o). 



(3.9) 



(3.10) 
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Figure 3. Solution of equation (3.121, with boundary condition v(.025) = 32.6 and parame- 
ter G — 67 ( ), corresponding to the values for figure[T]in dimensionless form. For comparison, 

the dashed line ( ) has G = 0. Both solutions become singular (v 1 —¥ — oo) at r s i ng = .85 

and 1.36, respectively. 



The two relations (3.9) can be used in (3.6) and in the mass conservation integral (3.5) 
to obtain 

vv' + ti = -G, vhr = l. (3.11) 
h 



For G = the equations (3.11) reduce to those of Bohr et al. (1993), which are essentially 
as derived by Kurihara ( 1946[ ) and Tani ( 1949 1. We combine (3.11 ) into one ODE for v(r): 

r(l - rv 3 )v' = -(1 - Gr 2 v - rV*)u , (3.12) 

which must be solved together with the flux boundary condtion v(ro) = «o &t the jet 
radius r . Note that there is a singularity at r — 0, and another at rv 3 = 1. The former is 
not relevant, since we have r > 0. The latter will determine the location of the hydraulic 
jump, which we will associate here with the rise distance. Since G — G(Q,8), this ODE 
will have to be solved at each inclination angle and flow rate. 

When doing numerical calculations, we will use the parabolic profile 

f(v) = h 2 -V 3 , (3.13) 
from which equation (3.10) gives c\ = 6/5, ci = 3. A more general approach, for instance 



using a variable cubic profile as in Bohr et al. (1997), doesn't significantly change the 
scaling. 



Let us examine solutions of (3.12) for typical experimental parameters. For the case 



shown in figure[T] Q ~ 119 cm 3 s and 9 = 45°. The viscosity of water is v ~ .01 cm 2 s~ 
and g ~ 980 cm s -2 . The jet radius is 0.28 cm, and from the flow rate this gives a 
velocity 484 cms -1 . Inserting all this into (3.4), we find a horizontal length scale T = 
11.3 cm, velocity scale a — 9.9 cms -1 , and G = 67. In dimensionless form, we must 



now integrate ( |3.12[ ) from r = (.28/11.3) = .025 with v Q = (2/3) (484/9.9) = 32.6. (The 
2/3 = 1//'(1) arises from the choice of a parabolic vertical profile for /.) 

Figure [3] shows the numerical solution (solid line), which becomes singular at r = 



.85. In dimensional form, the singularity is R s 



9.6 cm from the center of 



the jet. The experimentally-measured value for the lower edge of the hydraulic jump 



is i?;, 



7.2 cm, with is 3/4 of the singularity position. Note that these values are fairly 



insensitive to the exact jet radius r$. For comparison, the numerical solution for G = 
is also indicated as a dashed line in the figure, with a singularity at a larger value of r. 
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The singularity thus appears to occur somewhat beyond the actual position of the 
jump. To find the position of the jump, we would have to match to the 'outer' solution 



as done in Bohr et al 
parameter 



(small v) of (3.12), and impose continuity of mass and momentum across the jump, 

In addition, 



(1993) 



Bohr et al. 



( 1993 ) were left with an extra 
which they fixed by 



the location of a singularity of the outer solution 
assuming the jet was striking a plate of finite extent, and then moving the singularity 
to the edge of the plate. We cannot use this approach here: our hydraulic jump actually 
terminates, which by mass conservation means there must be either backflow or non- 
axisymmetric flow (or both). Both these effects require more powerful theories or the 
solution of more complex equations. 

On the other hand, since we are only after the scaling of the jump position, the theory 



we have is enough to uncover this scaling. Bohr et al. ( 1993 ) observed that their jump 
location typically occurred at unit radius (in dimensionless variables), somewhat inde- 
pendently of what was happening downstream in the outer solution. This suggests the 
following approach: fix the ratio of the jump distance to the singularity distance from 
the center of the jet. In our example above, that ratio was 3/4, but we find .76 fits the 



set of data slightly better. Hence, we have in dimensional form -Rj um p 
after using in T parabolic profile values for the numerical constant, c 



76r sing r, or 



1/2 -3/8 



.73: 



ump 



•55 r sins (G,r ,v ) x ^-^costf)" 1 / 8 . 



(3.14) 



Here r s i ng (G, ro, vq) is obtained by solving the dimensionless ODE (3.12| with initial 
condition at the jet radius v(ro) = Vq, and Vq is obtained from the dimensionless flow 
rate by vq = (2/3)(Q/7rrQ), where the 2/3 is for a parabolic vertical profile. 

In figure [2] we compare formula (3.14| with experiments at various inclination angles, 
with the jet always striking the plane normally. We emphasize that the numerical pref- 
actor in ( |3.14| is fixed, so we are not fitting each data set individually. The theory agrees 
well with experiments in both velocity and angle, except at low flow rates and at 90° 
angle. The low flow rate disagreement is not troubling, since it falls outside the theory as 



pointed out by Bohr et al. (1993). The 90° theoretical curve stands out, since it is very 



close to the theoretical curve for 45°. However, the trend of the 90° curve with Q is still 
captured. 

Figure [4] shows the singularity position as a function of angle, for a fixed flow rate 
119 cms -1 . The curve actually has a minimum around 80°, which is unphysical: physi- 
cally we expect the rise distance to decrease monotonically with angle, though at these 
large angles of inclination the flow is too unsteady for accurate measurement of the rise 
distance. We conclude that the theory breaks down for steeper angles: probably backflow 
becomes important, and the assumption that the the outer solution only weakly affects 
the jump position breaks down. According to figure |4j the rise distance ceases to increase 
significantly after about 60°; however figure [2] shows a significantly higher rise distance 
for 90° inclination. This suggests that the present theory works well for angles 9 < 60°. 



4. Conclusions 

We have presented simple experiments of the impact of a jet on an inclined plane, and 
noted the dependence of the maximum rise distance on both flow rate and angle. Though 
a complete description of this problem is daunting, our simple radial model captures the 
dependence remarkably well, though less so at smaller flow rates and larger angles. A 
complete theory would, of course, describe the thickness of the bounding rim, but also 



predict the critical angle at which the hydraulic jump changes from closed to open (Lebon 
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et aL| |2008). Our simple model verifies the need to include viscosity to capture the rise 
height at these modest flow rates, as pointed out in Edwards et al. (20081. 
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